function res = funcK(kh,nn1,nn2,dx,dy,dd)

ima = sqrt(-1);

res = ((nn1.*dx + nn2.*dy).*...
       (1./dd).*(-kh*besselh(1, kh*dd)) + ...
       ima*kh*besselh(0, kh*dd)...
       );

return